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ABSTRACT 

We present an optical spectroscopic analysis of the well-known Seyfert galaxy NGC 4051. The high- 
ionization nuclear emission-line region (HINER) traced by [Fe X] A6374 is found to be spatially extended 
to a radius of 3" (« 150 pc) west and southwest from the nucleus; NGC 4051 is the third example which 
has an extended HINER. 

The nuclear spectrum shows that the flux of [Fe X] A6374 is stronger than that of [Fe VII] A6087 in our 
observation. This property cannot be interpreted in terms of a simple one-zone photoionization model. 

In order to understand what happens in the nuclear region in NGC 4051, we investigate the physical 
condition of the nuclear emission-line region in detail using new photoionization models in which the 
following three emission- line components are taken into account; (1) optically thick, ionization-bounded 
clouds; (2) optically thin, matter-bounded clouds; and (3) a contamination component which emits Ha 
and H/3 lines. Here the observed extended HINER is considered to be associated with the low-density, 
matter-bounded clouds. Candidates of the contamination component are either the broad-line region 
(BLR) or nuclear star forming regions or both. The complexity of the excitation condition found in 
NGC 4051 can be consistently understood if we take account of these contamination components. 

Subject headings: galaxies: individual (NGC 4051) - galaxies: nuclei - galaxies: Seyfert 



1. INTRODUCTION 

It is known that Seyfert galaxies often show very high- 
ionization emission lines such as [Fe VII] A6087, [Fe X] 
A6374, [Fe XI] A7892 and [Fe XIV] A5303 (Oke & Sargent 
1968; Grandi 1978; Penston et al. 1984; De Robertis & Os- 
terbrock 1986). Because the ionization potentials of these 
lines are higher than 100 eV, much attention has been 
paid to the high-ionization nuclear emission-line region 
[HINER; Murayama, Taniguchi, & Iwasawa 1998 (here- 
after MTI98); see also Binette 1985]. The possible mech- 
anisms of radiating such high-ionization emission lines are 
the following three processes: (1) collisional ionization in 
the gas with temperatures of T^ ~ 10^ K (Oke & Sargent 
1968; Nussbaumer & Osterbrock 1970); (2) photoioniza- 
tion by the central nonthermal continuum emission [Os- 
terbrock 1969; Nussbaumer & Osterbrock 1970; Grandi 
1978; Korista & Ferland 1989; Ferguson, Korista, & Fer- 
land 1997b; Murayama & Taniguchi 1998a, 1998b (here- 
after MT98a and MT98b, respectively)]; and (3) a combi- 
nation of shocks and photoionization (Viegas-Aldrovandi 
& Contini 1989). 

Recently, in context of the locally-optimally emitting 
cloud models (LOG models; Ferguson et al. 1997a), Fer- 
guson, Korista, & Ferland (1997b) showed that the high- 
ionization emission lines can be radiated in conditions of 
widely range of gas densities. More recently MT98a have 
found that type 1 Seyfert nuclei (Sis) have excess [Fe VII] 
A6087 emission with respect to type 2s (S2s). Given the 
current unified model of AGN (Antonucci & Miller 1985; 
see for a review Antonucci 1993), the finding of MT98a 



implies that the HINER traced by the [Fe VII] A6087 emis- 
sion resides in the inner wall of such dusty tori. Since the 
covering factor of the torus is usually large (e.g., ~ 0.9), 
and the electron density in the tori (e.g., ~ lO"^"* cm~'^) 
is considered to be significantly higher than that (e.g., 
~ 1Q3-4 cm^"^) of the narrow- line region (NLR), the con- 
tribution from the torus dominates the emission of the 
higher-ionization lines (Pier & Voit 1995). Taking this 
HINER component into account, MT98b have constructed 
new dual-component (i.e., a typical NLR with a HINER 
torus) photoionization models and explained the observa- 
tions consistently. 

On the other hand, it is also known that some Seyfert 
nuclei have an extended HINER whose size amounts up 
to - 1 kpc (Golev et al. 1995; MTI98). The presence 
of such extended HINERs can be explained as the result 
of very low-density conditions in the interstellar medium 
(riH ~ 1 cni~^) makes it possible to achieve higher ioniza- 
tion conditions (Korista & Ferland 1989). Thus MT98a 
suggested a three-component model for the spatial distri- 
bution of HINER in terms of photoionization. That is: (1) 
the inner wall of the dusty torus with electron densities of 
no ~ lO^""^ cm-3; the torus HINER [Pier & Voit 1995; 
Murayama & Taniguchi 1998b], (2) the innermost part of 
the NLRs; the NLR HINER (n^ ~ lO^"" cm'^) at a dis- 
~^100 pc, and (3) the extended ionized 
cm~'^) at a distance ~1 kpc; the ex- 
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tended HINER (Korista & Ferland 1989; MTI98). Perhaps 
the relative contribution to the HINER emission from the 
above three components may be different from galaxy to 
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galaxy. In particular, extended HINERs have been found 
only in NGC 3516 (Golev et al. 1995) and Tololo 0109- 
383 (MTI 98) and thus it is important to investigate how 
common the extended HINER in Seyfert galaxies. 

In this paper, wc report on the discovery of an ex- 
tended HINER in the nearby Seyfert galaxy NGC 4051. 
This observation was made during the course of our long- 
slit optical spectroscopy program for a sample of nearby 
Seyfert galaxies at the Okayama Astrophysical Observa- 
tory. Throughout this paper, we use a distance toward 
NGC 4051 of 9.7 Mpc, which is estimated using a value of 
Ho = 75 km s^^ Mpc^^ and its recession velocity of 726 km 
s~^ (Ulvestad & Wilson 1984). Therefore, 1" corresponds 
to 47 pc at this distance. 

2. OBSERVATIONS 

The spectroscopic observations were made at Okayama 
Astrophysical Observatory, National Astronomical Obser- 
vatory of Japan on 1992 June 5. The New Cassegrain 
Spectrograph was attached to the Cassegrain focus of the 
188 cm reflector. A 512 x 512 CCD with pixel size of 24 x 
24 ^m was used, giving a spatial resolution of 1'.'46 pixel^^ 
by 1 X 2 binning. A I'.'S slit with a length of 300" was used 
with a grating of 150 groove mm~^ blazed at 5000 A. The 
position angle was set to 90°. The wavelength coverage 
was set to 4500 - 7000 A. We took three spectra; (1) the 
central region, (2) 2" north of the central region, and (3) 2" 
south of the central region. Each exposure time was 1200 
seconds, respectively. The slit positions for NGC 4051 are 
displayed in Figure 1. The data was reduced with the use 
of IRAF. The reduction was made with a standard proce- 
dure; bias subtraction and flat fielding were made with the 
data of the dome flats. The flux scale was calibrated by 
using a standard star (BD-|-33° 2642). The nuclear spec- 
trum was extracted with 2'.'92 aperture. The seeing size 
derived from the spatial profile of the standard star was 
about 2'.'3 (FWHM) during the observations. 

3. OBSERVATIONAL RESULTS 

3.1. Emission-Line Properties of the Nuclear Spectrum 

The spectrum of the nuclear region (the central 2'.'9 x 
1'.'8 region) is shown in Figure 2. In order to estimate 
emission-line fluxes, we made multicomponent Gaussian 
fitting for the spectrum using the SNG (the SpectroNeb- 
ularGraph; Kosugi et al. 1995) package. The identified 
emission lines of the nuclear region are summarized in Ta- 
ble 1. The [Fe X] A6374 emission line is blended with [O I] 
A6364. Assuming the theoretical ratio of [O I] A6300/[O I] 
A6364 = 3 (Osterbrock 1989), we measured the [Fe X] 
A6374 flux. The reddening was estimated by using the 
Balmer decrement (i.e., the ratio of narrow components of 
Ha and H/3). If the case B would be assumed, an intrinsic 
value of Ha/H/3 ratio was 2.87 for T = lO'' K (Osterbrock 
1989). However, Veilleux & Osterbrock (1987) mentioned 
that the harder photoionizing spectrum of AGNs results 
in a large transition zone or partly ionized region in which 
collisional excitation becomes important (Fcrland & Nct- 
zer 1983; Halpern and Steiner 1983). The main effect of 
the collisional excitation is to enhance Ha. Therefore we 
adopt Ha/H/3 ~ 3.1 for the intrinsic ratio, and accord- 
ingly we obtain Ay = 1.00 mag. This value is almost 
consistent with the previous estimation (Ay — 1.11 mag: 



Erkens et al. 1997). In our observation, [Fe X] A6374 (ion- 
ized potential 233.6 eV) is stronger than [Fe VII] A6087 
(99.1 eV). This observational result is inconsistent with 
the predictions of simple one-zone photoionization models 
(see section 4.). 

In Table 2, we give a comparison between our observa- 
tional data and the previous ones (Anderson 1970; Grandi 
1978; Yee 1980; Pension et al. 1984; Veilleux 1988; Erkens, 
Appenzeller, & Wagner 1997). Although Erkens et al. 
(1997) gave [Fe VII] A6087/[Fe X] A6375 = 0.966 in their 
paper, they newly reduced their data and found the true 
observed line ratio is 0.500 (Wagner & Appenzeller 1999, 
private communication). Though [Fe VII] A6087/[Fe X] 
A6374 in Veilleux (1988) is signiflcantly larger than that 
in ours, our ratio is consistent with those in Penston et al. 
(1984) and Erkens et al. (1997). Although we do not un- 
derstand fully the signiflcant difference between Veilleux 
(1988) and the other observations, it may be due partly 
to the difference of slit width or aperture size among the 
observations. 

NGC 4051 is one of the well-known Seyfert galaxies 
(Seyfert 1943). It has been mostly classified as a type 
1 Seyfert (Adams 1977), while Boiler, Brandt, & Fink 
(1996) and Komossa & Fink (1997) pointed out that the 
observational properties of NGC4051 are similar to those 
of narrow-line Seyfert 1 galaxies (NLSl; Osterbrock & Da- 
hari 1983; Osterbrock & Pogge 1985). Though our obser- 
vational data show the broad component of Ha clearly, 
the broad component of H/3 is not detected. The results 
of deconvolution for Ha and H/3 are shown in Figure 3 and 
4, respectively. 

3.2. Spatial Distribution of the Emission- Line Region 

In Table 3a - 3d, we give the emission-line properties 
of the off-nuclear regions; west (2'.'9 west), southeast (1'.'5 
south 2"west), southwest (1'.'5 south 2"east), and east (2'.'9 
east). Since the flux of [O I] A6300 in these areas is not 
measured because of the insufficient S/N, wc do not sub- 
tract the fiux of [O I] A6364 from that of [Fe X] A6374. 
Though wc measured the fluxes of emission lines of north- 
east, those data are not tabulated because we could not 
detect the H/3 unambiguously. The S/N of the northwest 
position is so poor that we did not measure the fluxes of 
emission lines. 

Figure 5 shows that the HINER traced by [Fe X] A6374 
is extended westward up to 3" (« 150 pc). This is more 
extended than the NLR traced by [O I] A6300. Since, as 
shown in Figure 6, there is no strong line of sky emission at 
the observed wavelength of [Fe X] , the extended [Fe X] ap- 
pears to be real. Figure 5 also shows that the HINER may 
be extended southwestward. However this may be due to 
the contamination from the nuclear region, suggested by 
the relatively broad width of H/3 at the southwest position. 

Following Veilleux & Osterbrock (1987), we investigate 
the excitation conditions of the emission-line region in 
each position. As shown in Figure 7, we flnd that the re- 
gions where [Fe X] is absent exhibit AGN-like excitations, 
whereas the regions where [Fe X] is found show H II region- 
like excitations (except for the southeast region where the 
line ratios show H II region-like excitation though [Fe X] is 
not detected) . It is unlikely that the [Fe X] emission arises 
from H II regions. Therefore the observed H II region-like 
excitations are due not to photoionization by massive stars 
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Fig. 1. — Slit positions wc set for NGC 4051. Images are taken from digitized sky survey. 
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Fig. 2.— Nuclear spectrum of NGC 4051. 
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Table 1 
Emission-Line Data of the Nuclear Region in NGC4051 



Identification 


'^obs 


FWHM=' 


F/F(RI3N)^ 




I/I{Ul3Ny 






(A) 


(km s-l) 


Line Ratios 


lad 


Line Ratios 


1 a" 


H/3N 


4874.6 


934.2 


1.000 


f 


1.000 




OIII 


4970.7 


431.3 


0.457 


±0.016 


0.445 


±0.015 


OIII 


5018.7 


427.1 


1.348 


±0.027 


1.294 


±0.026 


[Fe XIV] 






<0.07is 




<0.063S 




He I 


5886.5 


839.6 


0.166 


±0.014 


0.130 


±0.012 


Fe VII] 






<0.080S 




< 0.060S 




OI] 


6314.2 


453.7 


0.156 


±0.012 


0.113 


±0.009 


OI] 


6377.7 


449.2 


0.052 


±0.012 


0.037 


±0.009 


FeX] 


6387.4 


762.7 


0.190 


±0.014 


0.136 


±0.011 


Nil] 


6562.1 


458.9 


0.409 


±0.017 


0.286 


±0.014 


HaN 


6577.4 


666.3 


4.351 


±0.073 


3.030 


±0.075 


HqB 


6577.5 


3439.7 


2.285 


±0.040 


1.591 


±0.035 


[Nil] 


6597.4 


456.4 


1.208 


±0.029 


0.838 


±0.028 


[SII] 


6733.5 


h 


0.187 


±0.012 


0.127 


±0.009 


[SII] 


6747.9 


h 


0.201 


±0.012 


0.136 


±0.009 



^Corrected for the instrumental width. 

''_F{H/3N) = 1.460 X W'^^ erg s-^cm-^ 

"^The reddening-corrected relative intensities. We adopted Ay = 0.995. 
Accordingly, /(H/3N) = 4.269 X 10~13 erg s^^cm-^^ 

•* Estimated 1 cr-lcvel error for F/F(H/3N). 

''Estimated 1 cr-lcvcl error for ///(H/3N). 

^Estimated 1 tr-level error for _F{H/3N) is 2.451 X 10^-*^^ erg s~^cm~^. 

^3(7 upper limit. 

''Narrower than the measurable limit (instrumental width). 



Table 2 
A Comparison of Our Data with Previous Data 





Anderson 70 


Grandi 78 


Yoo 80 


PFBWW 84" 


Veillcux 88 


EAW 97*= 


Ours 


Date 


66~67? 


74? 


74 Apr. 21 


77 Apr. 4 


78 Mar. 25 
78 Apr. 8 


92 Jul. 


92 Jun.5 


Telescope 


Wilson 2.5m 


Lick 3m 


Hale 5m 


Hale 5m 


Lick 3m 


DSAZ'^2.2m 


GAG 1.8m 


Detector 


photomultiplicr 


ITS"^ 


photomultiplicr 


IPCS 


ITS'* 


CCD 


CCD 


Slit width 


slitless 


2". 7 


slitless? 


? 


2". 7 


1".5 


1".8 


[O IIIl/H/?" 


1.750 




2.399 








1.739 


[FoVII]/H/3f 




0.062 










<0.060 


[FoVII]/[FoX]s 








0.490 


1.000 
0.826 


0.500' 


<0.441 


[FoXI]/[FoX]'' 








0.324 




0.514' 





"Pcnston et al. (1984) 

''Erkens, Appenzcllcr, & Wagner (1997) 

^The Calar Alto Observatory (Spain) 

Image-Tube Scanner 
"[GlII] AA4959+5007/H/3narrow 
'[Fe VII] A6087/H/3narrow 
«[FeVII] A6087/[FeX] A6374 
•■[FeXI] A7892/[FeX] A6374 
'These value are obtained from EAW (1999; private communication). 
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Fig. 3. — Result of multi-Gaussian component fitting for the Ho-I- [N II] emission lines of NGC 4051. The wavelength is shown in unit of 
V, and the flux is shown in unit of erg cm~^ sec~^A. 
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Fig. 4. — Result of multi-Gaussian component fitting for the H/3+ [O III] emission lines of NGC 4051. The unit for each axis is the same 
as Figure 3. 
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Table 3a 
Emission-Line Data of the Off-Nuclear Region (2'.'9 west) 



Identification 



(A) 



FWHM'^ 
(km s^-*-) 



F/F(H/3N)'' 
Line Ratios 



///(H/3N)<: 
Line Ratios 



H/3N 
[O III] 
[O III] 

[FeX] 
[Nil] 
HaN 
Nil] 
SII] 
SII] 



4874.8 
4970.7 
5018.7 
6391.2 
6562.6 
6577.3 
6597.9 
6733.0 
6747.4 



677.1 
657.9 
651.6 
765.9 
771.2 
769.5 
767.1 

g 

g 



1.000 


t 


1.000 




0.563 


±0.068 


0.548 


±0.066 


1.660 


±0.115 


1.594 


±0.111 


0.680 


±0.075 


0.488 


±0.061 


0.634 


±0.075 


0.442 


±0.059 


4.440 


±0.271 


3.087 


±0.268 


1.870 


±0.128 


1.296 


±0.120 


0.268 


±0.053 


0.182 


±0.038 


0.405 


±0.056 


0.274 


±0.042 



''Corrected for the instrumental width. 

'^F{H/3N) = 2.216 X lO-l-* erg s'^cm'^ 

•^The reddening-corrected relative intensities. We adopted Ay = 1.000. 

•^Estimated 1 cr-lcvel error for observed line ratios. 
"Estimated 1 cr-level error for reddening-corrected line ratios. 
^Estimated 1 tr-level error for F(H/3N) is 1.315 X 10~^^ erg s~^cm~-^. 
s Narrower than the measurable limit (instrumental width). 



Table 3b 
Emission-Line Data of the Off-Nuclear Region (1'.'5 south - 2"west) 



Identification 



'^obs 

(A) 



FWHM'^ 
(km s^^) 



F/F(H/3N)b 
Line Ratios 



1 a^ 



///(H/3N)<^ 
Line Ratios 1 cr*^ 



H/3N 
[O III] 
[O III] 
He I 

[FeX] 
[Nil] 
HaN 

[Nil] 

[SII] 
[SII] 



4869.4 
4967.1 
5015.1 
5886.6 
6381.2 
6559.5 
6574.2 
6594.8 
6731.1 
6745.5 



839.7 
444.6 
440.3 
832.6 
705.4 
666.9 
665.4 
663.3 

g 

g 



1.000 


t 


1.000 




0.459 


±0.031 


0.443 


±0.030 


1.353 


±0.052 


1.285 


±0.049 


0.226 


±0.035 


0.166 


±0.026 


0.244 


±0.027 


0.159 


±0.018 


0.616 


±0.032 


0.388 


±0.024 


4.911 


±0.160 


3.083 


±0.143 


1.819 


±0.064 


1.137 


±0.055 


0.212 


±0.022 


0.129 


±0.014 


0.185 


±0.021 


0.112 


±0.014 



^Corrected for the instrumental width. 

''F(H/3N) = 4.378 x 10"" erg s^^cm^^ 

•^The reddening-corrected relative intensities. We adopted Ay = 1.280. 

'^Estimated 1 cr-level error for observed line ratios. 
''Estimated 1 cr-level error for reddening-corrected line ratios. 
^Estimated 1 cr-level error for F(H/3N) is 1.412 X 10"^^ erg s^-'cm"^. 
s Narrower than the measurable limit (instrumental width). 
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Table 3c 
Emission-Line Data of the Off-Nuclear Region (1'.'5 south - 2"east) 



Identification 



'^obs 

(A) 



FWHM'' 
(km s-^) 



Line Ratios 



1 a^ 



///(H/3N)<: 
Line Ratios 



1 cr-^ 



H/3N 
[O III] 
[O III] 

[Nil] 
HaN 

[Nil] 
[SII] 
[SII] 



4870.8 


1220.2 


4967.0 


589.4 


5015.0 


583.8 


6559.4 


674.7 


6574.2 


673.2 


6594.7 


671.1 


6729.7 


s 


6744.1 


R 



1.000 


t 


1.000 




0.406 


±0.056 


0.394 


±0.055 


1.196 


±0.091 


1.145 


±0.087 


0.627 


±0.054 


0.425 


±0.045 


4.566 


±0.290 


3.086 


±0.279 


1.849 


±0.122 


1.245 


±0.115 


0.245 


±0.027 


0.161 


±0.021 


0.287 


±0.029 


0.189 


±0.023 



^Corrected for the instrumental width. 

'>F(H/3N) = 3.492 x lO-l-* erg s-^cm-^ 

•^The reddening-corrected relative intensities. We adopted Ay = 1.078. 

"^Estimated 1 cr-level error for observed line ratios. 
''Estimated 1 cr-lcvcl error for reddening-corrected line ratios. 
^Estimated 1 cr-level error for F(H/3N) is 2.202 X 10^^^ erg s^^cm^^. 
s Narrower than the measurable limit (instrumental width). 



Table 3d 
Emission-Line Data of the Off-Nuclear Region (2'.'9 east) 



Identification 


'^obs 


FWHM^ 


F/F(H/3N)t' 




///(H/3N)<= 






(A) 


(km s-i) 


Line Ratios 


lad 


Line Ratios 


1 o-° 


H/3N 


4871.1 


781.7 


1.000 


f 


1.000 




OIII] 


4968.5 


517.5 


0.995 


±0.148 


0.979 


±0.146 


OIII] 


5016.5 


512.6 


2.936 


±0.338 


2.864 


±0.332 


Nil] 


6560.7 


607.2 


0.739 


±0.132 


0.594 


±0.126 


HaN 


6575.5 


605.8 


3.852 


±0.436 


3.092 


±0.498 


[Nil] 


6596.0 


604.0 


2.180 


±0.261 


1.747 


±0.291 



^Corrected for the instrumental width. 

^F(UI3N) = 1.720 X lO-l-* erg s-^cm-^ 

•^The reddening-corrected relative intensities. We adopted Ay = 0.604. 

"^Estimated 1 cr-level error for observed line ratios. 
''Estimated 1 cr-level error for reddening-corrected line ratios. 



^Estimated 1 cr-level error for F(H/3N) is 1.892 X 10 



-15 



erg s ^cm 
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but to some additional mechanism. We will discuss this 
complex property in section 5. 

3.3. A Summary of the Observational Results 

As noted in section 1, there are three kinds of the 
HINER; 1) the torus HINER, 2) the NLR HINER, and 3) 
the extended HINER (see MT98a). Our detection of the 
extended [Fe X] emitting region tells us that the extended 
HINER exists at least in NGC 4051. Here we estimate 
how strong the contribution from the torus HINER using 
the dual component photoionization model of MT98b. Ac- 
cording to the diagnostic diagram of their model (Figure 
2 in MT98b), we find that the torus HINER may con- 
tribute to the total intensity of the HINER emission less 
than 3%. On the other hand, if [Fe X] in the nuclear re- 
gion was mainly attributed to the NLR HINER, this line 
would have a larger FWHM in the nuclear region than 
in off-nuclear regions because the flux contribution of the 
NLR HINER may be negligibly small in the off-nuclear 
regions. Since we could not find the difference of FWHM 
of [Fe X] between the nuclear region and west, the NLR 
HINER is not a dominant source in NGC 4051. It is there- 
fore suggested that the majority of the HINER emission 
in NGC 4051 arises from low-density ISM within a radius 
of « 150 pc. 

4. PHOTOIONIZATION MODEL 

In order to understand the nuclear environment of NGC 
4051, we use photoionization models and compare the pre- 
dictions of the models with the observed emission-line ra- 
tios of the nuclear region of NGC 4051. The simplest 
model for the NLRs of Seyfert galaxies is a so-called one- 
zone model, which assumes optically thick clouds of single 
density and single distance from the source of the ion- 
ization radiation (e.g. Ferland & Netzer 1983, Stasinska 
1984). However, these models have been known to predict 
too weak high-ionization emission lines such as [Fe VII] 
A6087 and [Fe X] A6374, and moreover, predict less in- 
tense [Fe X] A6374 than [Fe VII] A6087. Because these 
model predictions appear inconsistent with observations, 
one-zone models are not suitable to investigate the envi- 
ronment of the nuclear region of NGC 4051. Hence, it 
is better to use more realistic models, for example, the 
optically-thin, multi-cloud model (Ferland & Osterbrock 
1986) or the LOG models (Ferguson et al. 1997a). The 
predicted emission-line flux calculated with these models 
are shown in Table 4. The former model predicts [Fe VII] 
A6087/[Fe X] A6374 « 10, that is inconsistent with the 
result of observations of NGC 4051. The LOG model 
also predicts [Fe X] < [Fe VII]. Taking these points into 
accounts, we construct a two-component photoionization 
model below. 

4.1. A Two- Component Photoionization Model 

We construct a two-component system for the nuclear 
emission-line region of NGC 4051. One component is op- 
tically thick, ionization-bounded clouds (IB clouds). This 
component emits low-ionization emission lines mainly, like 
clouds in typical NLRs. The other is optically thin, 
low-density, matter-bounded clouds (MB clouds; Viegas- 
Aldrovandi 1988; Viegas-Aldrovandi & Gruenwald 1988; 
Binette, Wilson, & Storchi-Bergmann 1996; Wilson, Bi- 
nette, & Storchi-Bergmann 1997; Binette et al. 1997) 



which radiates high-ionization emission lines selectively. 
These MB clouds are expected to emit more intense [Fe X] 
A6374 than [Fe VII] A6087 because their densities are as- 
sumed to be low enough to achieve very high-ionization 
conditions (section 4.2). 

Ionization and thermal equilibrium calculations have 
been performed with the photoionization code CLOUDY 
(version 90.04; Ferland 1996) to calculate the emission 
from plane-parallel, fixed hydrogen density clouds. Taking 
into account many lines of evidence in favor of a nitro- 
gen overabundance (Storchi-Bergmann & Pastoriza 1990; 
Storchi-Bergmann 1991; Storchi-Bergmann et al. 1998), 
we adopt twice the solar nitrogen abundance. Namely, all 
elements have solar values except for nitrogen. The detec- 
tion of strong [Fe X] suggests that most of iron remains 
in gas phase although the depletion of iron would be more 
serious than that of other elements (e.g., Phillips, Gond- 
halekar, & Pettini 1982). Therefore internal dust grains 
in the NLR are not taken into account in our calculations. 
The shape of the ionizing continuum from the central en- 



gine IS 



f„ = z/""" cxp(- 



hv 



BB 



)exp(- 



hv ' 



(1) 



We adopt the following parameters; (1) fcTiR is the infrared 
cutoff of the so-called big blue bump component and we 
adopt fcTiR = 0.01 Ryd; (2) Tbb is the temperature which 
parameterize the big blue bump continuum, and we adopt 
a typical value, 1.5 x 10^ K; (3) ctuv is the slope of the low 
energy big blue bump component. We adopt auv = —0.5. 
Note that the photoionization is not sensitive to this pa- 
rameter. And, (4) ttx is the slope of the X-Ray component, 
and we adopt Qx = —1.0. This power law component is 
not extrapolated below 1.36 eV or above 100 keV. Below 
1.36 eV, this term is set to zero while above 100 keV, the 
continuum is assumed to fall off as v~^ . Finally, (5) the 
UV to X-Ray spectral slope, aox, is defined as 



log[f^(2keV)/F,(2500A)] 
log[j/(2keV)/j/(2500A)] '' 



(2) 



which is a free parameter related to the parameter a in 
equation (1). We adopt aox — —1-4. The observational 
values of these parameters for NGC 4051 are summarized 
in Table 5. 

The calculations for IB clouds are proceeded until the 
electron temperature drop below 3000 K, since the gas 
with lower temperature than 3000 K is not thought to con- 
tribute significantly to the emission lines. The calculations 
for MB clouds are proceeded till the column density of the 
MB clouds reach to a value given as a free parameter. 

4.2. Results 

First, we discuss the physical conditions of the IB clouds. 
Assuming that the low-ionization forbidden lines are radi- 
ated mainly from the IB clouds, we estimate the hydrogen 
density of the IB cloud njB = 10^'^ cm~^, which derived 
from the observed [S II] doublet ratio, [S II] A6716/[S II] 
A6731 — 0.934 (see Osterbrock 1989). Similarly, assum- 
ing that the MB clouds contribute to the flux of the low- 
ionization lines very little, we search an ionization pa- 
rameter for the IB clouds, Um = (5(H)/(47r_R^A^H,iBc) 
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Table 4 
Model Predictions for Forbidden Lines 



Linc^ 


F086!' optically thin 


FKBF97^ LOG model 








(Solar Abundance) 


(Dusty Abundance)'* 


NeV A3426 


0.09 


0.62 


0.67 


O III A5007 


7.11 


11.7 


10.8 


Fe XIV] A5303 


C 


0.033 


0.004 


Fc VII] A6087 


0.02 


0.049 


0.009 


O I] A6300 


0.85 


0.56 


0.84 


Fc X] A6374 


0.002 


0.015 


0.002 


N II] A6583 


2.56 


1.18 


1.69 


S II] AA6716+6731 


1.37 


1.20 


1.87 


[Fe XI] A7892 


C 


0.068 


0.008 



^Normalized by H/3narrow. 

''Ferland & Osterbrock (1986) 

"^Ferguson et al. (1997a) 

•* Assuming the abundance of the Orion nebula (see FKBF97). 

''Not predicted by their model. 



Table 5 
Observational Constraint for Input Continuum 



parameter 



observational value adopted value 



Tbb 





1.5 X 10^ 


Q-uv 





-0.5=' 


Clx 


-1.88±0.18t' (ROSAT) 

-Om^iGinga) 

-0.85±0.07'= (ASCA) 

-1.30^1 (ROSAT) 


-1.0 


O-ox 


-1.32'' (ROSAT) 


-1.4<= 



''A recommended value in CLOUDY (See Ferland (1996) 
and Francis (1993).) 

''Walter et al. (1994) 

'^Guainazzi et al. (1996) 

"^Komossa & Fink (1997) 

''A recommended value in GLOUDY (See Ferland (1996) 
and Zamorani et al. (1981). ) 
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(the ratio of the ionizing photon density to the Hydro- 
gen density) using CLOUDY. Emission-hne ratios of [O I] 
A6300/[O III] A5007, [N II] A6583/[0 III] A5007 and [S II] 
AA6717,6731/[0 III] A5007 are calculated for various val- 
ues of UiB, and compared with the observed values. As 
shown in Figure 8, the comparisons for the individual line 
ratios do not give a certain value of Uib- Therefore, we use 
the [S II] AA6717,6731/[0 III] A5007 ratio to determine the 
hydrogen density of the IB clouds, and then we derive Uib 
_ xo~2-^ accordingly. 

Second, we estimate most probable values of the param- 
eters for MB clouds. When Umb > 10^"'^, the calculated 
flux of [Fe X] A6374 is smaller than that of [Fe XI] A7892; 
and when Umb < 10^°'^, [O III] A5007 begins to emit from 
MB clouds. Because these conditions are not suitable to 
explain the observations, we adopt Umb — 10"°-^. When 
the hydrogen column density of MB clouds A^mb > 10^^ 
cm^^, [O III] A5007 also begins to radiate from the MB 
clouds (see Figure 9). Therefore we examine two cases 
for A^MB — 10^°-^ cm^^ and for 10^^° cm^^. Assuming 
the size of HINER Z^hiner = 150pc == 4.63 x 10^° cm, 

10 



ionization emission line, we obtain another relation: 



we obtain umb — 10 "'^^ cm ^ for A^mb 



and Umb — lO"-^"' cm -^ for A^mb 



10" 



because 



riMB — A^mb/^'hiner- Since the former density is too low 
to produce sufficiently strong emission, we adopt the lat- 
ter case, that is, umb = 10°'^^ cm~'^ and A^mb = 10^^ 
cm^^. In Table 6, we give the emission line fluxes normal- 
ized by up (narrow component) for the IB and MB clouds 
described above. 

It seems reasonable that the nuclear emission-line region 
of NGC 4051 is a mixture of both IB and MB clouds. In 
order to reproduce the observed [Fe X]/H/3 ratio, we find 
that the relative contribution of the MB clouds is 5.3% in 
the H/3 luminosity. We compare the total calculated line 
ratios with the observed values in Table 7. We find that 
[O III] , [N II] and [S II] are two or three times stronger than 
the observational values although high-ionization lines are 
consistent with the observation. This discrepancy can be 
reconciled if there is another emission component which 
radiates hydrogen recombination lines mainly. Hereafter 
we call this "contamination component" . Possible con- 
tamination sources are either BLR or nuclear star forming 
regions or both (section 5). Note that this contamina- 
tion component must not destroy the line ratios such as 
AA6717,6731/[0 III] A5007, [S II] A6716/[S II] A6731 and 
so on. 

Flux contributions of the IB clouds, the MB clouds, and 
the contamination component to Balmer lines are treated 
as free parameters a and b; a is the H/3 flux ratio between 
the MB clouds and the IB clouds and b is that between the 
contamination component and the IB clouds. We can find 
a probable set of the line ratios which is consistent with 
the observation. Here we assume a ratio of Ha/H/3 for 
the contamination component to be 3.1. Because [Fe X] is 
assumed to emit from the MB clouds, we obtain a relation: 



, [Fe X ] , 
■ H/3 



Jobs 



--0^) 



MB 



1 + a + b 



(3) 



[o ni] 

H/3 



jobs 



( [O HI] 

^ H/3 



IB 



1 + a + b 



(4) 



Since we can regard [O III] A5007 as a representative low- 



In Table 6, we give ([Fe X]/H/3)ob«, ([Fe X]/H/3)mb, 
([O III]/H/3)obs, and ([O III]/H/3)ib. Using these relations, 
we find a = 0.161 and b = 1.868. These results mean that 
the contributions of the IB clouds, the MB clouds and the 
contamination component are 33.0%, 5.3% and 61.7% in 
the H/3 luminosity, respectively. We give a summary of the 
set of the line ratios in Table 8. Though we do not observe 
[Fe XI] A7892, previous observations (Pension et al. 1984 
and Erkens et al. 1997) show [Fe XI] A7892/[Fe X] A6374 
= 0.324 or 0.514. Our calculated [Fe XI] A7892/[Fe X] 
A6374 is 0.674, and this is not so different from the previ- 
ous observations. On the other hand, the observed [O I] 
A6300 is four times stronger than the model value. One 
reason for this may be that [O I] arises partly from other 
regions that we do not take into account in our model. 

Finally, in Table 9, we give a summary of the three 
emission components adopted for the nuclear emission-line 
region of NGC 4051. These parameters are determined 
uniquely in the process described above. However, there 
may be other models which explain the observed line ratios 
of NGC 4051. Recently, Contini & Viegas (1999) proposed 
a multi-cloud model in which the existence of shocks is in- 
troduced for NGC 4051. Their model explains the optical 
line ratios and the continuum SED, although they did not 
mention the spatial extension of ionized regions. In or- 
der to discriminate which model is more plausible, further 
detailed observations will be necessary. 

5. DISCUSSION 

As we have shown in previous sections, the observed 
emission- line ratios of the nuclear region of NGC4051 are 
consistently understood by introducing the three emission 
components; 1) the ionization-bounded clouds, 2) matter- 
bounded clouds, and 3) the contamination component to 
the Balmer emission lines. Although our three-component 
model appears consistent with the observation, our result 
implies that the majority of Balmer emission {'^ 60%) 
arises from the contamination component. Now we con- 
sider the problem; What is the contamination component 
? 

First, we consider this problem for the nuclear region. 
Possible candidates of the contamination components are 
either the BLR or nuclear star-forming regions or both. If 
NGC 4051 belongs to a class of NLSls (Boiler et al. 1996; 
Komossa & Fink 1997), it seems hard to measure the con- 
tribution from the BLR to the H/3 emission because of the 
narrow with of the broad line if present. It is also possible 
to consider that NGC 4051 experiences a burst of massive 
star formation in its nuclear region because there is a lot 
of cold molecular gas as well as circumnuclear star-forming 
regions in NGC 4051 (Kohno 1997; Vila-Vilaro, Taniguchi, 
& Nakai 1998). Peterson, Crenshaw, & Meyers (1985) re- 
ported that H/3 of NGC 4051 exhibits time variability (en- 
hanced by 85% in H/3 flux) on time scales shorter than ~2 
years, that is, H/3 of NGC4051 contains the broad compo- 
nent more or less. Although we have no way to evaluate 
the contribution of this BLR contamination to the total 
flux quantitatively, it is possible that all of the contamina- 
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Table 6 
Observed and Calculated Line Ratios 



Line 


Observed^ 


IB Clouds 


MB Clouds 


H/3N A4861 


1.000 


1.000 


1.000 


OIII A5007 


1.294 


3.937 


0.114 


Fc XIV] A5303 


< 0.063 


0.000 


0.160 


He I A5876 


0.130 


0.141 


0.001 


Fe VII] A6087 


< 0.060 


0.000 


0.014 


O I] A6300 


0.113 


0.091 


0.000 


Fc X] A6374 


0.136 


0.000 


2.658 


HaN A6563 


3.030 


2.901 


2.706 


N II] A6584 


0.838 


2.200 


0.000 


S II A6716 


0.127 


0.405 


0.000 


S II A6731 


0.136 


0.441 


0.000 


[Fc XI] A7892 




0.000 


1.784 



''The reddening-corrected relative intensities. 



/(H/3N) = 7.406 X 10~ 



erg s cm 



Table 7 
Comparison of Calculated Line Ratios'" IB+MB IVIodel with Observed Values 



Line 


IB clouds 


MB clouds 


Total 


Observed'' 


H/3N A4861 


0.947 


0.053 


1.000 


1.000 


O III] A5007 


3.728 


0.006 


3.734 


1.294 


Fc XIV] A5303 


0.000 


0.008 


0.008 


< 0.063 


He I A5876 


0.134 


0.000 


0.134 


0.130 


Fc VII] A6087 


0.000 


0.001 


0.001 


< 0.060 


O I] A6300 


0.086 


0.000 


0.086 


0.113 


Fc X] A6374 


0.000 


0.141 


0.141 


0.136 


HaN A6563 


2.747 


0.143 


2.890 


3.030 




Nil A6584 


2.083 


0.000 


2.083 


0.838 




S II A6716 


0.384 


0.000 


0.384 


0.127 




S II A6731 


0.418 


0.000 


0.418 


0.136 




Fc XI] A7892 


0.000 


0.095 


0.095 





^Normalized by observed H/3narrow. 

''The reddening-corrected relative intensities. 
/(H/3N) = 7.406 X lO^^^ erg s-^cm-^. 






[SII]X6716/[SII]A,6731 
^ (calculated value) 



[SII]A,6716/[SII]A,6731 
(observed value: 0.934) 




^ 



[Nil] >.6583 / [OIII]>. 5003 
(observed value: 0.648) 



[Sn]Mi6716+6731 / [OIII]>. 5003 
"■^^ (observed value: 0.203) 



[01] )l6300 / [Oin]X 5003 
(observed value: 0.087) 



log n(H) (cm ) 



logU 



Fig. 8. — Determination of the parameters for ionization-bounded clouds. The straight lines show the observed reddening-corrected line 
ratios and the curves show the calculated line ratio for various parameters (here, hydrogen density and ionization parameter). The small 
circles are located at the crossover for the lines describing the observed and calculated line ratios, therefore the circles indicate the estimated 
parameters. 
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Table 8 
Comparison of Calculated Line Ratios'* of Contaminated Model with Observed Values 



Line 


IB clouds 


MB clouds 


Contamination'' 


Total 


Observed*^ 


H/3N A4861 


0.330 


0.053 


0.617 


1.000 


1.000 


III] A5007 


1.299 


0.006 


0.000 


1.305 


1.294 


Fe XIV] A5303 


0.000 


0.008 


0.000 


0.008 


< 0.063 


He I A5876 


0.047 


0.000 


0.000 


0.047 


0.130 


Fc VII] A6087 


0.000 


0.001 


0.000 


0.001 


< 0.060 


OI A6300 


0.030 


0.000 


0.000 


0.030 


0.113 


Fe X] A6374 


0.000 


0.141 


0.000 


0.141 


0.136 


H«N A6563 


0.960 


0.143 


1.913 


3.016 


3.030 




N II] A6584 


0.726 


0.000 


0.000 


0.726 


0.838 




SII A6716 


0.134 


0.000 


0.000 


0.134 


0.127 




S 11 A6731 


0.146 


0.000 


0.000 


0.146 


0.136 




Fe XI] A7892 


0.000 


0.095 


0.000 


0.095 





'^Normalized by observed H/3narrow. 

''We assume an ratio of Ha/H/3 for this contamination component to be 3.1. 

"^The reddening-corrected relative intensities. 
/(H/3N) = 7.406 X 10"'^ erg s^'cm-^. 



Table 9 
Properties of Each Component 



(cm^^) 



U Nii Contribution to H/3 

(cm-2) 



IB clouds 
MB clouds 
Contamination 



102.9 
100.33 



10-2.9 
10-0.5 



10^ 



33.0% 
5.3% 
61.7% 



''The hydrogen density of the IB clouds is not a free parameter in our 
calculation. This value is actually enough high to make IB clouds optically- 
thick. 



log U = -0.6 



log U = -0.5 



log U = -0.4 








^ "^^^-^--^ [OIII] _ 


[Fe X] / 


I^\^,^^ 


[Fe XI[ ~Z 


y 


' ^--^ ^"^"---^ 


"[Fe^vll] 




[FeXIV] 



22.5 20.5 




log N(H) (cm"^ ) 



Fig. 9. — Determination of the parameters for matter-bounded clouds. The curves show the calculated fluxes normalized by H/3 for the 
varying hydrogen column density in the case of log V = —0.6, —0.5 and -0.4. 
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tion component is contributed from the BLR. In addition, 
the nuclear star-formation may contribute to the contami- 
nation. Kohno (1997) discussed the gravitational instabil- 
ity of the nuclear molecular gas of some Seyfert galaxies 
using the Toomre's Q-value. The Toomrc's Q parameter 
characterizes the criterion for local stability in thin isother- 
mal disks and is expressed as Q ^ Scrit/Sgas, where Scrit 
is the critical surface density. He gave Q — 0.90 for the 
nuclear region of NGC4051. This means that the molecu- 
lar gas in the nuclear region of NGC4051 is thought to be 
gravitationally unstable. 

In any case, about 60% of observed H/3 is not originated 
from the NLR in the nuclear region of NGC 4051. This 
means that the line ratios of the nuclear region suffer seri- 
ously from the contamination. In Figure 10, we replot the 
excitation diagnostic diagram using the line ratios, from 
which the contamination component is subtracted. This 
diagram shows that the contamination-subtracted line ra- 
tios of nuclear region show the typical AGN-like excita- 
tion condition. Therefore we conclude that the unusual 
excitation condition is due to the contamination compo- 
nent. High-spatial resolution optical spectroscopy or X- 
ray imaging observations will be helpful in investigating 
whether or not the star-formation activity dominates the 
flux of H/3. 

Second, we consider the off-nuclear regions. As shown 
in Figure 7, the three off- nuclear regions (west, southwest, 
and southeast) also show H II region-like excitations. Since 
a typical size of the BLR is ~ 0.01 pc (e.g. Peterson 1997), 
it is likely that these excitation conditions are thought to 
be due to circumnuclear star- forming regions. 
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log([NII] A6583)/(HaX6563) 



Fig. 10. — A diagnostic diagram saincly as Figure 6. The filled circle shows the excitation condition of the nuclear region corrected for the 
contamination effect. 
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